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Abstract. Modeling type la supernova (SN la) explosions in three dimensions allows to eliminate 
any undetermined parameters and provides predictive power to simulations. This is necessary to 
improve the understanding of the explosion mechanism and to settle the question of the applicability 
of SNe la in cosmological distance measurements. Since the models contain no tunable parameters, 
it is also possible to directly assess their validity on the basis of a comparison with observations. 
Here, we describe the modeling of SNe la as thermonuclear explosions in which the flame after 
ignition near the center of the progenitor white dwarf star propagates outward in the sub-sonic 
deflagration mode accelerated by the interaction with turbulence. We explore the capabilities of this 
model by comparison with observations and show in a preliminary approach, how such a model can 
be applied to study the origin of the diversity of SNe la. 
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INTRODUCTION 

Type la supernovae (SNe la) have become one of the major tools in observational cos- 
mology. Yet a sound theoretical understanding of these objects - justifying in particular 
the calibration techniques applied in distance measurements - is still lacking. First at- 
tempts to model SNe la were based on one-dimensional numerical simulations. Such 
models gave valuable insight into the basic mechanism of SN la explosions. However, 
their predictive power is limited due to the fact that underlying physical processes en- 
ter the models in a parametrized way. This is overcome by three-dimensional modeling 
of SNe la JHI2]. As an example, we present a model that is derived from the standard 
scenario of SN la explosions (for a review see Jj]). A white dwarf (WD) consisting of 
carbon and oxygen is assumed to accrete matter from a non-degenerate binary compan- 
ion until its mass approaches the Chandrasekhar limit. Due to the rapid increase of the 
central density nuclear reactions ignite giving rise to a stage of convective carbon burn- 
ing. This stage lasts for several hundred years and terminates once the nuclear energy 
production cannot be balanced by convective cooling any longer. Subsequently, a ther- 
monuclear runaway of a small temperature fluctuation ignites a thermonuclear flame. 
The exact mechanism of flame ignition, however, remains controversial. While some 
studies suggest a flame ignition in multiple sparks distributed around the center of the 
WD 0,BlIa0], others put forward central single-point ignitions 10]. 

After ignition the flame propagation is determined by the laws of hydrodynamics. 
Regarding the flame front as a discontinuity between fuel and ashes, they allow for two 
distinct modes of flame propagation. In a sub-sonic deflagration burning is mediated 
by thermal conduction while a super- sonic detonation is driven by a shock wave. On 



the basis of one-dimensional simulations of a prompt detonation, Arnett [9] ruled out 
this model for SNe la since it drastically underproduces intermediate mass elements 
observed in spectra of these events. 

Starting out as a laminar deflagration, however, the flame propagates too slowly to ex- 
plain the energy release necessary to explode the WD. Thus, any valid SN la model needs 
to provide means of flame acceleration. Two mechanisms are conceivable here. Firstly, 
the flame propagation may continue in the deflagration mode being significantly accel- 
erated by the interaction with turbulence. The one-dimensional model W7 of Nomoto 
et al. [10] demonstrated that such a model is in principle capable of reproducing the 
main observational features of SNe la. An alternative way to speed up the flame is to as- 
sume a deflagration-to-detonation transition (DDT) at later stages of the explosion. The 
weak point of these delayed detonation models (e.g. WllIP is that a physical mechanism 
providing a DDT in SNe la could not be identified yet lfl2l 12, 14, 3] and therefore 
the hypothetical transition of the burning mode enters the model as an undetermined 
parameter. 



A DEFLAGRATION TYPE la SUPERNOVAE MODEL 

Our goal in the following is to present a SN la explosion model that contains no tunable 
parameters. Therefore we set aside the possibility of a delayed detonation and focus on 
the turbulent deflagration model. Turbulence is induced here by generic instabilities. The 
flame propagates from the center of the star outward leaving behind light and hot ashes. 
Dense and cold fuel in front of the flame gives rise to an inverse density stratification 
in the gravitational field of the WD. This renders the flame propagation buoyancy 
unstable and in its nonlinear stage the Rayleigh-Taylor instability leads to the formation 
of (typically mushroom-shaped) burning bubbles that rise into the cold fuel. At the 
interfaces of these bubbles strong shear flows emerge. The corresponding Reynolds 
numbers reach values of the order of 10 14 and therefore strong turbulence is generated 
by secondary shear (Kelvin-Helmholtz) instabilities. The turbulent eddies generated 
on scales of the buoyancy-induced flame features decay to smaller scales forming a 
turbulent energy cascade and interact with the flame propagation. This stretches and 
corrugates the flame enlarging its surface area and thus the net burning rate is increased. 

The main challenge in numerically implementing this scenario is the vast range of 
relevant length scales involved. Not only is the width of a thermonuclear flame in the 
degenerate carbon/oxygen material at the onset of the explosion 9 orders of magnitude 
below the radius of the WD. The turbulent cascade extends to even smaller scales 
and interacts with the flame down to the Gibson length at which the laminar flame 
speed equals the turbulent velocity fluctuations (10 4 cm and decreasing in the explosion 
process). This problem can be tackled in a Large Eddy Simulation (LES) approach. Here, 
only the largest scales of the problem are directly resolved (applying the Prometheus 
implementation lllfll for solving the hydrodynamics equations)and turbulence effects 
on unresolved scales are included via a subgrid-scale model Ill7tll8ll . The thermonuclear 
flame is modeled as a sharp discontinuity separating the burnt from the unburnt material. 
Its evolution is followed utilizing the level-set method IU9I1 . Since the structure of the 
flame is not resolved, the flame propagation velocity must be provided externally. This, 




FIGURE 1. Snapshots from a full-star SN la simulation starting from a multi-spot ignition scenario. 
The logarithm of the density is volume rendered indicating the extend of the WD star and the isosurface 
corresponds to the thermonuclear flame. The last snapshot marks the end of the simulation and is not on 
scale with the earlier snapshots. 



however, does not introduce an undetermined parameter to the model since the theory 
of turbulent combustion [20] predicts that for most stages of the SN la explosion the 
flame propagation proceeds in the flamelet regime where it completely decouples from 
the microphysics of the burning and is determined by the turbulent velocity fluctuations 
that can be derived from the subgrid-scale model. This implementation provides a self- 
consistent model of SN la explosions in the deflagration scenario. 

Our description of burning is augmented by a simplified treatment of the nuclear 
reactions including only five species lElll . It provides the energy release necessary to 
follow the explosion dynamics. In order to derive observables from the models, however, 
the chemical composition of the ejecta needs to be known in detail. This is achieved in a 
nucleosynthesis postprocessing step 12211 . Of particular interest is the yield of 56 Ni, since 
its radioactive decay powers the light curve. 



TABLE 1. Variation of initial parameters in SN la explosion models 



Parameter 


Range of variation 


Effect on 56 Ni production 


Effect on total energy 


X( 12 C) 


[0.30, 0.62] 


<2% 


-14% 


p c [10 9 g/cm 3 ] 


[1.0,2.6] 


-6% 


-17% 


Z[Z Q ] 


[0.5, 3.0] 


-20% 


none 



SIMULATION RESULTS 

A typical evolution of a SN la explosion modeled as described above is shown in 
Fig. [T] Starting from an ignition in multiple sparks the flame propagates outward. At 
t = 0.3 s, the mushroom-shaped features due to the buoyancy instability are clearly 
visible. Subsequently, the flame becomes increasingly corrugated and is accelerated 
by interaction with turbulence. It therefore burns through a large fraction of the WD 
material. The snapshot at t = 0.6 s shows the flame evolution around the peak of energy 
production due to nuclear burning. Up to this point, the burning terminated in nuclear 
statistical equilibrium (NSE) and the carbon/oxygen material was primarily converted to 
iron group elements. The expansion of the WD decreases the fuel density steadily and 
once it falls below 5 x 10 7 gcm -3 nuclear burning becomes incomplete and produces 
mainly intermediate mass elements. About 2 s after ignition, expansion quenches the 
burning and the following evolution is characterized by the relaxation to homologous 
expansion of the ejecta, which is reached to a reasonable accuracy ~10s after ignition 
||23J]. The density structure of the ejecta at this stage is shown in Fig.d where the traces 
of turbulent flame propagation are clearly visible. 

Apart from the initial conditions simulations as described above contain no free pa- 
rameters. Therefore the question arises whether such models are capable of reproducing 
observations without any fitting. The explosion energies achievable in the outlined sce- 
nario reach up to ~8 x 10 50 erg and the models produce ~0.4 M Q of 56 Ni. This falls into 
the range of observational expectations, although on the side of the weaker SN la explo- 



sions [24, 25]. Nonetheless, the synthetic lightcurves derived from models of the class 
described here fit the observations in the B and V bands around maximum luminosity 
rather well A much harder constraint on the explosion model is posed by spec- 

tral observations, since spectra are particularly sensitive to the chemical composition of 
the ejecta. Kozma et al. 112 811 pointed out a potential problem of deflagration SN la mod- 
els. In late time "nebular" spectra, unburnt material (transported towards the center in 
downdrafts due to the large-scale buoyancy-unstable flame pattern) gives rise to a strong 
oxygen line of low-velocity material which is in conflict with observations. However, 
the synthetic spectrum of [28] was derived from a simplistic centrally ignited model. Re- 
cently, Ropke et al. l29ll showed that multi-spot ignition models may succeed to burn out 
the central parts of the WD, reducing the amount of oxygen at low velocities. A stochas- 
tic multi-spot ignition leads to similar results [30]. Detailed spectral observations allow 
to determine the chemical composition of the ejecta in velocity space ll3lll . The mixed 
composition of the ejecta observed there points to a deflagration phase being at least a 
significant contribution to the SN la explosion process. The central parts are found to 
be clearly dominated by iron group elements, which are mixed out to velocities of about 




1.4 1.5 1.6 1.7 1.4 1.5 1.6 1.7 

Am 15 (5) Ami 5 (fl) 



FIGURE 2. Peak luminosity vs. decline rate of the light curve in the B band (diamonds correspond to 
SN la explosion models; the dashed curve in the left panel indicates the original relation by Phillips et al. 
1 34] and a shifted relation is marked as a solid curve) 



1 2000km s . Intermediate mass elements are distributed over a wide range in radii and 
no unburnt material is found at velocities <5000kms _1 . A recent high-resolution full- 
star deflagration SN la simulation Il32ll demonstrated that such models can get close to 
these observational constraints. 

Although it cannot be ruled out that the pure deflagration model of SNe la is in- 
complete, the results indicate that it may be at least a dominant part of the mechanism. 
Therefore it is justified to ask how such models are affected by initial parameters of the 
exploding WD. This may give a hint to the origin of the diversity of SNe la. To mod- 
erate the computational expenses, simplified setups may be used to study the effects of 
physical parameters on the explosion models. Such an approach was recently taken by 
Ropke et al. [33] and resulted in the first systematic study of progenitor parameters in 
three-dimensional models. The basis of this study was a single-octant setup with moder- 
ate (yet numerically converged) resolution. However, the lack in resolution did not allow 
for a reasonable multi-spot ignition scenario and thus only weak explosions can be ex- 
pected. It is therefore not possible to set the absolute scale of effects in this approach, 
but trends can clearly be identified. The parameters chosen for the study were the WD's 
carbon-to-oxygen ratio, its central density at ignition and its 22 Ne mass fraction resulting 
from the metallicity of the progenitor. All parameters were varied independently to study 
the individual effects on the explosion process. The results of this survey are given in 
Tab.[T] To determine the effects of these variations on observables, synthetic light curves 
were derived from all models 1220 . From these, the peak luminosities and decline rates 
(in magnitudes 15 days after maximum; Am^) were determined. A comparison with the 
relation given by Phillips et al. l34ll (forming the basis of the calibration of cosmological 
distance measurements) is provided in the left panel of Fig. El Obviously, the absolute 
magnitude of the Phillips relation is not met by our set of models. Moreover, the range 
of scatter in Am^ is much narrower than that of the set of observations used by Phillips 
et al. H. But there is a trend of our models consistent with the slope of the Phillips 
relation. The right panel of Fig. El shows that this slope is dominated by the variation of 



the progenitor's metallicity. 
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